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Abstract. We present the first optical polarimetric observations of the runaway microquasar LS 5039. Our results reveal the 
presence of a large amount (~5%) of polarized emission towards this binary system. By combining polarimetric and spectro- 
scopic observations of some stars in the field together with available statistical information on the galactic interstellar polar- 
ization of the region, we have estimated and subtracted the contribution of the interstellar polarization in this direction. As a 
result, we obtain an intrinsic polarization of ~3% for the object, much higher than what would be expected from jet emission in 
the optical domain. We suggest that the polarized light originates by electron Thomson scattering in the stellar envelope of the 
companion star. This allows us to constrain the size of the region where the polarized emission originates, as well as estimating 
the scattering electronic density and the wind velocity at such distance. 
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1. Introduction 

The high mass X-ray binary system LS 5039 is one of the about 
sixteen confirmed Galactic microquasars (Paredes et al. 2002 
Rib6'5n02i. Recent astrometric studies carried out by Ribd 
et al. (2002 1, show that it is a runaway system escaping from 
the Galactic plane with a systemic velocity of ~150 km s ' . 

LS 5039 is a bright y~11.2 star with an ON6.5V((f)) 
spectral type (McSwain et al. ,2004), located at a distance of 
2.9 + 0.3 kpc (Ribo et al. l2002> . The binary system has a short 
orbital period of Po^b = 4.4267 + 0.0005 days and a high ec- 
centricity of e = 0.48 + 0.06 (McSwain et al. 2004 1. At ra- 
dio frequencies, the source has a non-thermal spectral index 
with moderate, but not periodic, variability. No strong radio 
outbursts have been ever detected (Marti et al. .l998. Ribd et al. 

rr^ Ribo Euna . 

In X-rays, the source presents a hard spectrum up to 30 keV, 
but no pulsed nor periodic emission has been detected (Ribo 
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et al. Reig et al. UM^ . The possibihty that L S 5039 

is a y-ray emitter was suggested by Paredes et al. <2000> . 
who proposed its association with the unidentified y-ray source 
3EG J1824-1514 (Hartman et al. [19991. In this context, the y- 
rays would be produced by inverse Compton upscattering of 
lower-energy photons by the non-thermal relativistic popu- 
lation responsible of the radio emission. Based on wind accre- 
tion models McSwain et al. (2004i suggest a neutron star as 
the compact object in the binary system, which would have an 
inclination angle between 40° and 60°. 

Although several physical parameters of the system are 
known, there is not yet a clear picture allowing to explain the 
broadband emission detected from it. Therefore, new obser- 
vations that help to constrain the properties of the system are 
needed. 

Optical polarimetric observations are a useful tool in the 
study of intrinsic polarization of X-ray binaries, since they can 
give us important information about their physical and geo- 
metrical properties. This kind of studies have been carried out 
by Scaltriti et al. <1997> for the microquasar GRO J 1655 -40, 
where they found a significant amount of intrinsic linear polar- 
ization (~3%) in the VRI bands, with the polarization direction 
being parallel to the accretion disk plane. In a following paper 
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on this object, Gliozzi et al. (I1998> detected oscillations in the 
polarization consistent with the orbital period. 

The intrinsic polarization in X-ray binary systems can re- 
sult from nonthermal emission processes due to the presence of 
high-energy electrons in the relativistic jets or, alternatively, by 
electron Thomson scattering due to an extended plasma above 
an accretion disk or a stellar envelope. Moreover, this intrinsic 
polarization could change with the binary phase, due to the po- 
larizing mechanism, to the distribution and physical state of the 
material, or by geometric factors such as the inclination angle. 

The main problem in the effective identification of intrinsic 
polarization of a source at optical frequencies is the contam- 
ination produced by the Galactic interstellar medium (ISM) 
polarization. This component is thought to arise due to ex- 
tinction by anisotropic dust particles preferentially oriented 
by the Galactic magnetic field. At present, the most complete 
study of Galactic starlight polarization has been carried out by 
Fosalba et al. f'2002V They presented a statistical analysis of 
the Galactic ISM polarization using ~5500 stars and found that 
the interstellar grains are not fully aligned with the Galactic 
magnetic field, a fact that was interpreted as the effect of a 
large random component of the field. However, it is important 
to stress that an entire reconstruction of the three-dimensional 
magnetic orientations requires additional radio, submillimeter 
and infrared data. 

In this paper, we report the first optical polarimetry of 
LS 5039. We describe the observations in Sect. |2] present the 
obtained results in Sect. |3l and we discuss on the origin of the 
detected polarization and state our conclusions in Sect.|4] 

2. Optical observations 

2.1. Polarimetry 

Polarimetric observations of LS 5039 were obtained at the 
2.15 m "Jorge Sahade" telescope, CASLEO, Argentina, us- 
ing a two-channel polarimeter This instrument is a rotating 
plate polarimeter with a Wollaston prism that splits the inci- 
dent light beam into two components, each one directed to a 
different photomultiplier (Martinez et al.p^990l. Its optical de- 
sign is based on VATPOL, the Vatican Observatory polarimeter 
(Magalhaes et al..l984j. 

The observations spanned four consecutive nights, roughly 
centered at the new moon, from 2003 June 28 to July 1 (UT 
dates). Weather conditions were fairly good for most of the ob- 
serving time, except for thick cirrus during the second night 
and thin cirrus during a few hours at the beginning of the first 
and third nights. The target was repeatedly observed through 
V and / filters, closely matching the Johnson-Cousins system, 
with individual exposure times of 300-360 s; this yielded be- 
tween 8 and 12 individual data points each night in each band. 
Sky integrations were regularly obtained between science ex- 
posures. A 11. 3" diaphragm was used for all the observations. 

Six field stars, with apparent magnitudes bracketing that of 
LS 5039 and within ~3.5' from its position (see Fig.0, were 
also observed in order to estimate the foreground polarization. 
Both polarized and unpolarized stars from the lists of Turnshek 
et al. J1990t were also observed each night, in order to evaluate 
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Fig.l. DSS 7"x7" image centered on LS 5039, showing the 
target and the six comparison stars. Polarization vectors in the 
V band have been overlaid. North is up and East is to the left. 



the instrumental polarization (which was found to be practi- 
cally zero) and to correct the measured polarization angles to 
the equatorial standard system. We were able to determine this 
zero point correction to an accuracy of +5°. 

Since both the ordinary and the extraordinary rays are mea- 
sured almost simultaneously, polarimetric observations should 
be barely affected by non-photometric conditions. However, 
large and/or rapid variations in the sky polarization (e.g., when 
the Moon is rising or setting) are expected to introduce system- 
atic errors in the measured polarization, because the sky condi- 
tions may have changed between the object and sky measure- 
ments. To prevent against this source of error, we interpolated 
the sky brightness and polarization (in the Stokes U - Q plane) 
to the same instant of each object observation. 

As a further check, we plotted the magnitudes and Stokes 
parameters of both object and sky as a function of time, in or- 
der to verify whether or not any change in the observed target 
polarization may have been caused by a sudden change in the 
sky polarization and/or transparency. This allowed us to reject 
some suspicious data points, as will be explained below. 

2.2. Spectroscopy 

Spectroscopic observations for the most polarized (~5%) com- 
parison stars were conducted on 2003 October 13 with the 
2.56 m Nordic Optical Telescope (NOT) in the island of La 
Palma (Spain). The instrument used was the Andalucia Faint 
Object Spectrograph and Camera (ALFOSC) combined with 
Grism 6 and the CCD detector EEV42-40 (2Kx2K). The CCD 
frames were processed using standard tasks in the IRAF pack- 
age, including bias subtraction, flat-fielding, optimal extraction 
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Table 1. Average polarization and position angle measure- 
ments of the comparison stars and LS 5039. 



Object Name 


(Pv) 




(Pi) 




USNO-A2.0 0750- 


(%) 


n 


(%) 


n 


13157416 (#1) 


6.36±0.41 


155.2+1.8 


4.75+0.11 


155.9+0.7 


13153701 (#2) 


0.32±0.04 


67.7±3.7 


0.42+0.01 


58.5+1.0 


13157934 (#3) 


4.92±1.00 


137.5+5.8 


5.67+0.08 


162.3+0.4 


13162766 (#4) 


0.58±0.02 


54.2+1.2 


0.45+0.02 


48.7+1.3 


13170022 (#5) 


0.51 ±0.04 


65.2+2.0 


0.47+0.03 


57.2+1.7 


13155723 (#6) 


0.47+0.07 


54.8+4.3 


0.40+0.06 


58.1+3.9 


LS 5039 


5.26±0.15 


161.3+0.8 


4.70+0.16 


162.4+1.0 



and wavelength calibration based on arc frames taken with a 
helium lamp. 

3. Results 

The polarization and position angle measurements in the V and 
/ bands of the six comparison stars, and average values for 
LS 5039, are quoted in TableQ] We show in Fig.^a Digitized 
Sky Survey (DSS) image of the target field with polarization 
vectors in the V band overlaid onto LS 5039 and the six com- 
parison stars. 

It is evident that, at least, two different polarization patterns 
are produced by interstellar dust along the studied direction. 
Stars #2, #4, #5, and #6 show polarization vectors with small 
amplitudes ((f i/>~0.3-0.6%), roughly ahgned in the NE-SW 
direction ({9v}~6Q°). On the other hand, light from stars #1 
and #3 is polarized almost at an orthogonal direction {{0v)~ 
135-155°), with a significantly higher percentage ({Pv)~5- 
6%); hence, these two stars have polarization properties very 
similar to that of LS 5039. 

We show in Fig. 12 the polarization measurements of 
LS 5039 in the V (filled circles) and / (open circles) 
bands as function of the orbital phase (fo=HJD 2452756.49, 
forb=4.4267 d). During the first night, corresponding to 
<D-0.05 in Fig. |2j the sky polarization remained fairly con- 
stant at ~1.5%, both in V and /. Since, in addition, its flux 
contribution within the diaphragm was ~4.5-5 mag lower than 
that from LS 5039, no significant systematic errors from sky 
polarization are expected. The / band polarization shows no 
significant variations with time, remaining nearly constant at 
{Pi) - 4.6%. In turn, the V band polarization time-curve shows 
fluctuations at a ^3cr level, except for two points falling signif- 
icantly below the mean. These data correspond to brightenings 
in the light-curve, and the larger variation in Py was accom- 
panied by a ~4° change in 6y. Although instrumental errors as 
the sources for these variations cannot be absolutely ruled out, 
there is no indication that any such error could have induced 
the observed behaviour. However, this would contrast with the 
apparent photometric stability of LS 5039 within few hun- 
dredths of magnitude reported by several authors (e.g. Marti 
et al. l2004l . 
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Fig. 2. Polarization and position angle in the V (filled circles) 
and / (open circles) bands as a function of the orbital phase for 
LS 5039. 



On the second night (<1)^0.30), as stated above, we ob- 
served through thick cirrus. As a consequence, the quality of 
these data is not as good, with large fluctuations (up to ~4%) 
in the sky polarization; hence, significant systematic errors 
could be aff'ecting this data set. The first half of the third night 
(0^0.50) was affected by thin cirrus, so we had to discard a few 
suspicious data points. Finally, on the fourth night (<I)^0.75) we 
observed mostly under photometric conditions. The object was 
constantly ~4 mag brighter than the sky within the diaphragm, 
and no rapid changes in the sky polarization were evident. No 
significant systematic errors should therefore be expected. 

Concerning the NOT spectroscopic observations, the sky 
conditions were non-photometric and affected by clouds. 
Moreover, we could only acquire spectra with the objects be- 
ing at a significant air mass. Consequently, a spectrum suitable 
for classification was only obtained for the brightest compari- 
son star #3. After combining two 600 s exposures, the resulting 
spectrum is consistent with an AO lb star taking into account 
the strong narrow hydrogen absorption lines. 

On the other hand, the apparent magnitudes of star #3 are 
mB=13.92+0.04 and my=12.86±0.04, hence B - y=1.06+0.06 
(Marti et al. 120041 . For an AO supergiant, the absolute mag- 
nitude is My=-5. 3+0.1 and the intrinsic color {B - y)o=0.00 
(Allen rT973> . Therefore, this comparison star has a color ex- 
cess of £(5- y)=1.06±0. 06 magnitudes. Using the relationship 
Ay = (3.30 + 0.28(B - V)q + 0.04£b_v)£(B - V) (Schmidt- 
Kaler Q5H2ll' this corresponds to an interstellar extinction of 
Av=3.5+0.2 and a distance of 8.4+0.1 kpc. 
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4. Discussion 

The observed polarization in astrophysical objects is usually 
the sum of two components, the intrinsic one produced in the 
object and the one arising from the foreground ISM between 
the object and the observer The high polarization measured 
towards LS 5039 (~5% in V band) suggests that, at least, part 
of it should have an intrinsic origin because it is very unusual 
to be produced in the ISM up to a distance of 3 kpc. 

To deduce the intrinsic polarization of the source. Pin, it is 
necessary to know the contribution of the ISM polarization, i.e. 
the value of Pism and subtract it from the observed one. 

First of all, we assume that light emitted from stars #1 to 
#6 is not polarized. The similar percentages and angles of po- 
larization of stars #2, #4, #5 and #6, suggest that they are being 
slightly affected by the ISM and that they lie close to the Sun, 
possibly in the local spiral arm (distance <500 pc). On the con- 
trary, stars #1 and 3# have polarization properties similar to that 
ofLS 5039. 

As we have mentioned in previous sections, the main prob- 
lem to determine the intrinsic polarization of an object, is the 
absence of a complete study on the Galactic interstellar polar- 
ization. However, we can obtain a rather crude estimation of 
the intrinsic polarization of LS 5039 by combining the statis- 
tical study of galactic interstellar polarization carried out by 
Fosalba et al. ('2002') and all the information gathered by us 
from polarimetric and spectroscopic observations of the other 
stars in the field. 

Fosalba et al. ("20021, have found that the distribution of 
the polarization degree as a function of distance grows linearly 
up to ~2 kpc, but beyond this distance the behavior of P(%) is 
best fitted by third-order polynomials up to a distance of ~6 kpc 
(see their Fig. 3). Adopting a distance of ~2.9 kpc for LS 5039 
(Ribo et al. I2002t the polarization percentage as a function of 
the distance results ~2.3%. From the equations that give the 
behavior of the polarization degree with Galactic longitude and 
latitude, we obtain values for P(%)~0.8 and P(%)~2, respec- 
tively. On the other hand, if we assume that the contribution of 
the interstellar polarization increases linearly with distance, we 
derive Pism~2.5% a value consistent with those obtained from 
Fosalba et al. (2002'), for the set of known stellar parameters 
of LS 5039. So, we suggest that the intrinsic polarization of 
LS 5039 is -3%. 

This value is also similar to the one obtained by Gliozzi 
et al. ( ,1998) for the microquasar GRO J 1655 -40, which was 
interpreted as the consequence of electron scattering by plasma 
above the accretion disk. 

The origin of the intrinsic polarized optical emission es- 
timated for LS 5039 could be due to a nonthermal emission 
process or to electron Thomson scattering. If the observed po- 
larized flux is the result of synchrotron emission originated 
from the relativistic plasma inside the jet, then the radio elec- 
tron populations should play an important role in the gener- 
ation of optical emission. However, we have extrapolated the 
energy distribution of the jet in LS 5039 (Ribo 2002| from ra- 
dio to optical frequencies and find that it is at least 4 orders of 
magnitude lower than the observed flux. Since the nonthermal 
emission process is not enough to generate the polarized opti- 



cal emission detected towards LS 5039, the electron Thomson 
scattering appears as a suitable alternative mechanism causing 
the observed polarized flux. Since there is no clear evidence 
of the presence of a large accretion disk in LS 5039, the most 
likely possibility is that the polarized emission originates in the 
stellar envelope of the companion star Such scenario has been 
previously proposed by Brown & McLean (1977) to explain 
the intrinsic polarization by Thomson scattering on free elec- 
trons in optically thin stellar envelopes. An improved treatment 
is given in Brown et al. ( 1978), while other works on stellar 
scattering polarization can be found in Cassinelli et al. ( I1987> 
and references therein. 

Considering the Brown & McLean picture for the gener- 
ation of intrinsic polarization, we can estimate some physical 
and geometrical parameters for LS 5039. Assuming that the 
plasma fills an ellipsoidal shell of any thickness with a uniform 
density distribution no of scattering electrons, the expected po- 
larization is given by: 

P = o-7-noP(l - 3r)sin^! (1) 

where, ctt is the Thomson cross-section, R is the radius of the 
envelope, / is the inclination angle of the equatorial plane and 
y is a function of the ratio a of the equatorial to polar radius. 
Assuming a simple radial flow, the corresponding electron den- 
sity of the envelope is: 



47rP2v^[l-^J /zmn 

Following McSwain et al. ( I2004t . suitable values for LS 5039 
are a wind velocity at infinity Vco=2440+190 km s ', y6=1.0, 
an optical star radius Popt=8.5 Rq and a wind mass loss 
AJ^opt=4xlO"'* Mq yr"' (the last two quantities are mean val- 
ues). The other symbols have their usual meanings. Introducing 
Eq. (1) into Eq. (2) and considering values of y of 0.1, 0.2 
and 0.3 (P(%)>0) we estimate that the density of the scatter- 
ing electrons is in the range (0.2-1.3)xl0'2 cm"^. From the 
above results we conclude that the polarized emission orig- 
inates mainly close to the star at a distance between (8.51- 
8.6) Rq, where the velocity of the stellar wind is (3-20) km s"' . 

An interesting feature in the polarization of LS 5039 is the 
detected variability during the first and second nights. If the 
polarimetric fluctuations observed on timescales of hours are 
real, then they might be attributed to localized density enhance- 
ments (or blobs) propagating in the general wind. The existence 
of these blobs has been introduced by Brown et al. ( 1995^ to 
explain transient features seen in polarimetric light curves of 
WR stars. Further simultaneous photometric and polarimetric 
observations of the system are necessary to shed light on the 
origin of the variable optical polarization. 
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